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1.0  INTRODUCTION 


This  report  describes  progress  on  a  research  program  on  the  fundamental  processes 
in  partially  ionized  plasmas  conducted  in  the  High  Temperature  Gasdynamics  Laboratory  at 
Stanford  University.  This  research  is  supported  by  a  grant  from  the  Air  Force  Office  of 
Scientific  Research  (AFOSR-88-0264)  and  is  currently  conducted  under  the  direction  of 
Professor  Charles  H.  Kruger.  Two  Ph.D.  candidates  are  currently  finishing  their  doctoral 
research  and  three  students  have  completed  their  Ph.D.  under  this  program.  These  three 
currently  have  faculty  positions  at  other  universities.  Another  student  has  also  accepted  a 
faculty  position. 

Several  space  power  and  propulsion  systems  of  potential  long  range  interest  to  the 
Air  Force  involve  partially  ionized  plasmas.  Such  systems  include  MPD  thrusters,  both 
open  and  closed  cycle  MHD  power  generation,  high  velocity  reentry  and  thermionic  energy 
conversion.  Although  the  specific  configurations,  the  exact  operating  conditions,  and 
which  of  the  competing  systems  will  prove  to  be  most  useful  in  the  long  term  remain  to  be 
established,  it  is  important  at  this  time  to  provide  a  broad  fundamental  research  base  in 
support  of  development  activity.  In  particular,  there  are  a  number  of  key  issues  regarding 
the  properties  and  discharge  behavior  of  partially  ionized  plasmas  and  the  interaction  of 
discharges  with  fluid  dynamics  that  need  to  be  understood  before  the  potential  and 
limitations  of  competing  systems  can  be  fully  evaluated.  In  addition,  it  is  important  that 
outstanding  young  applied  scientists  be  educated  in  these  areas. 

Our  research  on  partially  ionized  plasmas  was  initiated  under  granc  AFOSR-83- 
0108  and  focused  on  three  major  areas: 

1 .  plasma  properties, 

2.  discharge  effects:  plasma  electrode  interaction, 

3 .  interaction  of  discharges  and  fluid  dynamics. 

Extensive  reports  on  areas  (2)  and  (3)  were  given  in  ‘ne  previous  Annual  Scientific 
Reports.  Research  this  year  has  emphasized  area  (1),  does  this  report. 

Progress  during  the  past  year  on  plasma  properties  is  described  in  the  following 
two  sections.  Publications  and  presentations  resulting  from  this  work  are  cited  in  Section  4 
and  Section  5  lists  the  personnel  who  have  contributed  to  this  report. 
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2-0  MEASUREMENTS  OF  TEIE  VOLUMETRIC  RADIATIVE  SOURCE  STRENGTH  OF 

AN  AIR  PLASMA  BETWEEN  5000  AND  7500  K 


1.  Summary 

The  volumetric  radiative  source  strength  of  an  atmospheric  pressure  air  plasma  has  been 
measured  between  5000  and  7500  K  in  a  50  kW  inductively  coupled  plasma  torch.  The  radiation 
source  strength  and  the  temperature  profiles,  across  the  axisymmetric  5  cm  diameter  plasma, 
have  been  determined  independently  and  combined.  Local  thermodynamic  equilibrium  (LTE) 
temperatures  (based  on  the  absolute  intensity  of  five  atomic  oxygen  and  nitrogen  lines), 
Boltzmann  temperatures  (based  on  the  relative  intensity  of  two  oxygen  lines),  and  rotational 
temperatures  from  the  relative  intensity  of  N^  rotational  lines,  have  been  measured.  All  these 
temperatiu-es  agree  within  experimental  uncenainty.  In  addition,  electron  number  densities, 
determined  from  the  width  of  the  Hp  line,  are  close  to  equilibrium  predictions.  Since  the 
foregoing  measurements  support  the  assumption  of  approximate  LTE  in  the  plasma,  the  radiation 
source  strength  measurements  may  be  compared  with  theoretical  predictions.  Good  agreement  is 
obtained  with  the  latest  version  of  the  NEQAIR  code,  but  earlier  models  overestimate  the 
radiation  by  up  to  a  factor  three.  The  discrepancies  may  be  attributed  to  inaccurate  spectroscopic 
data  and  partition  functions  in  the  different  models. 


2.  Nomenclature 

Einstein  emission  coefficient  (s'^). 
gj  Degeneracy  of  the  excited  level. 

Iji  Total  line  intensity  (W/cm^/sr). 

nj  Number  density  of  j*  level. 

P  Pressure. 

Qei  Electronic  partition  function. 

Sr  Volumetric  radiative  source  strength  (W/cm^sr). 
Tr  Boltzmann  temperature. 

Tue  Local  thermodynamic  equilibrium  temperature. 
Xs  Mole  fraction  of  species  s  (oxygen  or  nitrogen), 
ejj  Energy  difference  between  the  levels  i  and  j. 

X  Wavelength. 
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To  understand  and  model  high  temperature  flows  for  re-entry,  plasma  processing  or  other 
applications,  accurate  estimates  of  the  volumetric  radiative  source  strength,  or  local  emissivity, 
are  of  importance.  In  the  present  paper,  we  repon  measurements  of  several  plasma  parameters 
including  LTE,  Boltzmann  and  rotational  temperatures,  electron  number  density  and  radiation 
source  strength  between  5000  and  7500  K.  As  will  be  discussed  in  Section  5,  the  temperature 
and  electron  number  density  measurements  show  that  the  plasma  studied  here  is  close  to  LTE. 

We  are  only  aware  of  source  strength  measurements  above  9000  K.  These  were  made 
either  in  shock  tubes  (Page*,  Nerem^)  or  in  stabilized  arc  discharges  (Schreiber,^  Ogurtsova,^ 
Devoto^).  The  emitted  radiation  was  recorded  by  means  of  thermopiles  and  calibrated  with 
tungsten  ribbon  lamps.  Temperatures  were  determined  under  the  assumption  of  LTE  firom  the 
absolute  intensity  of  atomic  oxygen  and  nitrogen  lines.  The  spectral  range  covered  generally 
extended  from  0.2  |im  to  several  pm.  In  addition,  Schreiber^  extended  his  measurements  to  the 
UV  out  to  0.05  |im  using  an  energy  balance  method  incorporating  a  Hall  probe. 

In  the  early  1960s,  space  programs  prompted  the  development  of  theoretical  models  for  high- 
temperature  air  radiation.^'^  More  recently.  Space  Shuttle  re-entry  modeling  and  the  Mars 
Mission  sparked  new  interest  in  such  calculations.^  While  Kivel,^  Meyerott,^  and  Breene^'s 
models  are  for  air  plasmas  in  local  thermodynamic  equilibrium,  the  NEQAIR  code,  developped 

by  Park,9  can  also  account  for  thermal 
non-equilibrium  by  incorporating 
translational,  rotational,  vibrational  and 
electronic  temperatures.  For  air  in  LTE, 
significant  differences  exist  between  the 
results  of  these  four  models.  These 
discrepancies  are  attributable  to 
differences  in  the  methods  of  calculation 
and  in  the  radiative  transition 
4000  6000  8000  10000  12000  14000  probabilities.  A  comparison  with  the 

Temperature  [K]  present  radiation  source  strength 

measurements  is  presented  in  Section  6. 

Figure  1 .  Air  source  strength 


The  experiments  presented  here  were  conducted  with  a  nominally  50  kW  TAFA  model  66 
RF  induction  plasma  torch,  powered  by  a  LEPEL  model  T-50  power  supply.  The  vertical  5-coil 


torch  has  a  7.5  cm  inner  diameter  and  an 
inner  length  of  26.5  cm  measured 
between  the  gas  injection  plates  and  the 
5  cm  diameter  exit  nozzle,  which  is  6  cm 
above  the  upper  coil.  All  measurements 
were  made  1  cm  downstream  of  the 
nozzle  exit. 

Plasma  line  intensities  (Fig.2)  were 
measured  using  a  SPEX  model  1400-11 
3/4  meter  scanning  monochromator 
fitted  with  a  Hamamatsu  model  R1104 


Hg.2.  Experimental  sehip  lor  emissiori  measurements  photomultiplier  tube.  Absolute  intensity 

calibrations  were  obtained  by  means  of  a 

tungsten  strip  lamp,  with  a  calibration  traceable  to  NBS  standards.  The  brighmess  temperature 
was  measured  with  an  optical  pyrometer.  The  optical  arrangement  permitted  lateral  scanning  of 


the  plasma.  Data  were  acquired  using  a 
Stanford  Research  Systems  model 
SR510  lock-in  amplifier,  then  transferred 
and  stored  on  a  laboratory  computer  for 
processing. 

Total  radiation  measurements  (Fig.  3) 
were  made  using  a  self-calibrating 
Scientech  model  36-0001  disc 
calorimeter,  and  a  mirrorless  traversing 


system  including  a  UV  grade  quartz  lens. 


Fig.  3.  Experimental  set-up  for  radiation  source  strength 


5.  Measurements 


5.1.  Temperature  measurements 

For  a  plasma  at  atmospheric  pressure  the  assumption  of  local  thermodynamic  equilibrium  is 
often  assumed  to  be  satistied.  However,  for  purposes  of  verification,  LTE,  Boltzmann  and 
rotational  temperatures  have  been  determined.  In  addition,  we  also  measured  electron  number 
densities. 


5.  La.  LTE  temperature 

We  measured  the  total  intensity  from  atomic  oxygen  lines  at  777.3  and  615.7  nm,  and  from  a 
triplet  of  atomic  nitrogen  lines  at  742.3, 744.2  and  746.8  nm.  After  subtraction  of  the  continuum 
signal  and  calibration,  the  data  were  Abel-inverted.  The  absolute  line  intensity  temperatures, 
Tlte*  were  calculated  from  the  Maxwell-  Boltzmann  relation 

ni  XsP  1  z'  Ej 

where  nj  is  related  to  the  measured  line  emission  coefficient  Ijj  by 


^  ^  4tc  Xji 
j  he  Aji  ’ 


(2) 


and  where  Xs  stands  for  the  mole  fraction  of  O  or  N  in  air  at  chemical  equilibrium.  The  quantity 
XgPAtT  has  been  substituted  in  Eq.  (1)  for  the  atomic  density.  An  iterative  calculation  eventually 
provides  Tlte.  The  uncertainty  is  typically  2%  over  most  of  the  plasma  extent.  At  the  edges  it 
increases  to  about  6%  for  the  777. 3nm  O  line  (most  intense  line)  and  10%  for  other  lines. 


5.1.b.  Boltzmann  temperature 

For  lines  of  the  same  atom  with  sufficient  upper  state  energy  separation  (2eV  for  the  two 
aforementioned  oxygen  lines),  one  can  infer  a  temperatme  independently  of  the  ground  state 
population,  considering  only  that  the  populations  of  the  emitting  levels  are  in  a  Boltzmann  ratio: 


[lii  =  111 

8j2 


exp 


Combining  Eqs.  (2)  and  (3): 

Jjiii  ^  gji  ^ji»i  ^j2i2 


I 


J2»2  gj2  ^j2i2  ^jl«l 


ex 


■  ^12) 
kTB  ; 

(3) 

(4) 

M  kTs  / 

from  which  Te  is  readily  obtained.  Since  the  energy  separation  between  the  two  levels  is 
relatively  small,  the  results  are  quite  sensitive  to  the  uncertainty  on  experimental  data  and  on  the 
Einstein  spontaneous  emission  coefficients.  The  uncertainty  in  the  Boltzmann  temperatures  is 
typically  10%,  except  at  outer  radii  where  it  is  higher. 
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5.I.C.  Rotational  temperature 

The  rotational  lines  of  the  B  X  (0,0)  band  of  were  resolved.  They  yielded,  after 
Abel  inversion,  a  radial  rotational  temperature  distribution.  ^0  The  accuracy  of  the  resulting 
rotational  temperature  profile  is  however  rather  poor  since  the  uncertainty  on  lateral 
measurements  is  greatly  amplified  by  the  Abel  inversion.  For  this  reason,  only  a  temperature 
inferred  from  the  rotational  profile  measured  along  a  diameter  of  the  plasma  has  been  indicated 
in  Fig.4.  This  "average"  temperature  is  consistent  with  LTE. 


5.1.d.  Interpretation  of  the  temperature  results  (see  Figure  4) 

The  values  of  Tlte  different  lines  at  a  fixed  radial  position  are  typically  within  lOOK. 
The  Boltzmann  temperature  T3  agrees  within  experimental  uncertainty,  and  the  averaged 

rotational  temperature  is  comparable 
with  Tlte  and  supports  the 

assumption  of  approximate  LTE  of  the 
plasma  in  the  considered  temperature 
range.  Moreover,  since  the  atomic 
oxygen  and  nitrogen  line  intensities 
depend  (although  weakly)  on  the  degree 
of  dissociation,  the  data  are  consistent 
with  chemical  equilibrium  as  well. 

0  0.5  1.0  1.5  2.0  2.5 


Radial  position  [cm] 


5.2.  Electron  number  densit 


Figure  4.  Temperature  profiles 
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o  6.0x10'" 


4.0  X  10' 


2.0  X  10' 


Plate  Power  =  67  kW 
Air  Flow  Rate  =  94.9  l/min- 
nX  Hj  Flow  Rate  =  1  l/min 


—  Measured  Electron  Density 
—  Equilibrium  Electron  density 


Radial  position  [cm] 


Electron  number  densities  have  been 
measured  from  the  width  of  Abel- 
inverted  profiles  of  the  Hp  line.  About 
1%  of  H2  was  mixed  with  air  for  that 
purpose.  The  LTE  temperature 
measured  from  the  oxygen  line  at  777.3 
nm  dropped  uniformly  over  the  plasma 
extent  by  about  200  K.  Based  on  this 
temperature,  LTE  electron  number 
densities  were  calculated  with  the  Saha 
equation.  The  good  agreement  between 
these  calculations  and  the  measurements 


Figure  5.  Electron  density  profile 
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(Fig  5)  adds  further  support  to  the  LTE  assumption  over  most  of  the  plasma. 

5.3.  Radiation  source  stremth  measurements 

The  radiation  source  strength  as  a  function  of  lateral  position  was  measured  as  described  in 
Section  4.  The  data  were  Abel-inverted  and  reduced  to  volumetric  radiative  source  strength. 

Figure  6  includes  a  plot  of  the  radiation  source  strength  versus  the  LTE  temperature  obtained 
from  the  most  intense  line  (01 777.3  nm),  since  it  allows  for  the  best  accuracy. 


6.  Comparison  with  theoretical  predicrions 


As  mentioned  earlier,  calculations 
of  the  radiation  source  strength  have 
been  performed  by  several  authors 
under  the  assumption  of  LTE.  The 
predictions  of  Kivel^  and  Meyerott^  are 
higher  than  our  measurements  by  a 
factor  of  two  to  three.  This  discrepancy 
is  likely  caused  by  the  inaccurate 
radiative  transition  probabilities  then 
available.  Kivel  and  Meyerott's 
transition  probabilities  are  an  order  of 
magnitude  greater  than  those  used  by 
Breene.8  Breene's  predictions  are  much 
closer  to  our  measurements,  but  the  agreement  might  be  fortuitious  for  anomalies  in  the 
calculations  have  been  reported  by  Keck,ll  and  for  the  radiative  transition  probabilities  are 
somewhat  different  from  the  currently  accepted  values. 

Modifications  have  been  made  to  the  NEQAIR  code^  in  the  form  of  updated  radiative 
transition  probabilities,  and  of  partition  functions  recalculated  following  Lavrov. 1“*  The 
resulting  predictions  are  in  exccellent  agreement  with  our  measurements.  It  should  be  noted  that 
the  radiation  from  CN  has  not  been  included  in  the  NEQAIR  calculations  but  preliminary 
estimates  indicate  that  it  might  increase  the  total  radiation  by  up  to  10%  at  6(X)0  K. 


Figure  6.  Air  source  strength 
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7.  Conclusion 


Source  strength  measurements  have  been  conducted  in  a  plasma  found  to  be  close  to  LTE.  This 
plasma  appears  as  an  excellent  test  ground  for  high-temperature  air  radiation  modeling.  While 
early  models  tended  to  overestimate  the  radiation,  the  latest  version  of  the  NEQAIR  code 
reproduces  accurately  the  present  measurements  over  a  broad  temperature  range.  Spectral 
measurements  are  currently  under  way  for  purposes  of  comparison  with  NEQAIR  spectral 
calculations. 
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3-0  ELECTRONIC  QUENCHING  OF  ARGON  EXClltU  STATES  IN  A  NON-EQUILIBRIUM 

PLASMA  AT  ATMOSPHERIC  PRESSURE 


Summary 

Emission  measurements  of  temperature  and  electron  density  have  been  made  downstream  of  a 
50  kW  induction  plasma  torch  at  temperatures  and  electron  densities  near  7000  K  and  10^0  m*^, 
respectively.  Results  indicate  that  for  pure  argon  plasmas  the  free  and  bound-excited  electrons  are  in 
mutual  equilibrium  although  over-populated  with  respect  to  the  ground  state.  The  addition  of  either  a 
hvdroeen  or  nitrogen  diluent,  however,  disturbs  this  mumal  equilibrium  because  of  an  exothermic 
electronic  interaction  between  the  diluent's  atomic  ground  state  and  an  argon  excited  state. 

2.  Introduction 

Most  plasma  emission  diagnostics  rely  on  an  assumption  of  either  local  thermodynamic 
euuiiibnum  (LTE)  or  partial  LTE.  For  an  LTE  plasma,  only  one  temperature,  along  with  the 
pressure,  is  required  to  specify  the  thermodynamic  state,  while  for  a  PLTE  plasma  another  parameter 
;s  necessary.  A  PLTE  plasma  is  characterized  by  having  the  bound-excited  and  free  electrons  in 
mutual  equilibrium  though  over  or  under  populated  with  respect  to  the  ground  state.  The  amount  of 
this  over  or  under  population  defines  the  additional  parameter  required  to  specify  the  state  in  a  PLTE 
plasma.  Experimentally,  it  is  often  found  that  although  LTE  does  not  exist  in  many  atmospheric 
argon  plasmas.  PLTE  does  hold,/l-2/ 

While  pure  argon  plasmas  provide  an  excellent  and  relatively  simple  test  condition  for  plasma 
liaHnostics,  plasma  mixtures  provide  a  more  important  condition  enabling  the  study  of  plasma 
jhemistrv.  .Additionally,  mixtures  are  often  used  for  diagnostic  measurements  of  rotational 
•.emperature  from  molecular  transitions  or  of  electron  densities  from  Stark  broadened  transitions. 
Accordingly,  there  is  a  need  to  examine  the  applicability  of  various  emission  diagnostics  to  plasma 
mixtures. 

In  this  paper,  we  report  absolute  and  relative  atomic  line  intensities,  absolute  recombination 
continuum  in  both  the  visible  and  the  ultraviolet.  Stark  broadening  of  the  hydrogen  Balmer  lines,  and 
independent  calorimetric  measurements  which  demonstrate  that  the  presence  of  either  nitrogen  or 
hvdrogen  in  amounts  as  small  as  0.1%  cause  depanures  from  PLTE  in  a  recombining  atmospheric 
argon  plasma. 

V  Experimental  Facility 

The  measurements  desenbed  here  were  conducted  with  a  nominally  50  kW  TAFA  model  66  RF 
inducuon  plasma  torch,  powered  by  a  LEPEL  model  T-50  power  supply.  The  venical  5-coil  torch 
lias  a  7.5  cm  inner  diameter  and  an  inner  length  of  26.5  cm  measured  between  the  gas-injection  plate 
and  the  5  cm  exit  nozzle.  The  induction  coil  has  a  9  cm  diameter  and  a  7.6  cm  length.  The  present 
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expenmenis  were  conducted  at  two  locations.  Spectroscopic  measurements  reponed  at  the  "nozzle 
exit"  were  made  1  cm  downstream  of  the  acmal  5  cm  nozzle  exit  which  is  itself  7  cm  downstream  of 
the  last  induction  coil.  Spectroscopic  measurements  reponed  at  the  "test  section  exit"  were  made  1 
cm  downstream  of  a  water  cooled  quanz  test  section  which  has  an  inner  diameter  of  5  cm  and  a 
leneih  of  17.5  cm  above  the  nozzle  exit.  The  power  supply  and  torch,  and  the  test  section  cooling 
water  systems,  were  separately  instrumented  with  thermocouples  and  flowmeters  to  independently 
obtain  calorimetric  energy  balances  for  comparison  with  those  obtained  from  spectroscopic 
measurements. 

Plasma  emission  measurements  were  made  using  a  SPEX  model  1400-11  3/4  meter  scanning 
monochromator  fitted  with  a  Hamamatsu  model  R1104  photomultiplier  tube.  Absolute  intensity 
calibrations  in  the  visible  were  obtained  by  means  of  a  tungsten  strip  lamp,  with  a  calibration 
•raceable  to  NBS  standards.  Temperature  measurement  of  the  tungsten  strip  lamp  was  made  using  a 
['’yro-micro  disappeanng  filament  optical  pyrometer.  Ultraviolet  calibration  was  accomplished  using 
an  ’equilibrium"  condition  of  the  plasma  torch./3/  A  four  mirror  optical  arrangement  permitted  lateral 
and  axiai  scanning  of  the  plasma.  Data  sets  acquired  using  a  Stanford  Research  Systems  model 
SR51()  lock-in  amplifier  were  transferred  and  stored  on  a  laboratory  computer  for  processing.  The 
lateral  traverses  of  emission  were  transformed  to  radial  variations  within  the  plasma  using  an  Abel 
inversion  technique. 

4,  Experimental  Methods 


If.  as  asserted,  neither  LTE  nor  PLTE  prevails  in  our  plasma  mixtures,  the  experimental 
diacnostics  which  are  implemented  must  directly  measure  desired  parameters  independent  of  these 
equilibrium  assumptions,  .\bsolute  atomic  line  emission  measurements  lead  directly  to  the 
i^opuiaucn  m  m  ^  panicular  snecies  m  s  excited  state  through  the  following  relation: 


where  L,  is  the  emission  coefficient,  Ajj  is  the  transition's  wavelength,  h  is  Planck's  constant,  c  is 
the  speed  of  light,  and  Ajj  is  the  Einstein  transition  probability.  This  measurement  is  independent  of 
LTE  or  PLTE.  requiring  only  that  the  emission  is  optically  thin. 

.\bsolute  emission  measurements  of  recombination  continuum  allow  one  to  evaluate  both 
absolute  electron  number  densities  as  well  as  electron  kinetic  temperatures  which  are  also  both 
indepenaent  of  LTE  or  PLTE./4/  In  general,  the  recombination  emission  depends  on  the  electron 
number  density,  n^,  the  ion  number  density,  oj.  the  electron  temperature.  Te,  the  wavelength,  X.  and 
the  lias  composition.  Radiative  recombination,  the  inverse  process  of  photoionizaiion.  is 


symbolically  represented  as  follows: 


.VI-*-  -i-  e  M  +  hv 
hv  =  ejx  -I-  Ext 


(2) 
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The  energy  of  the  emitted  photon,  hv,  is  the  sum  of  the  energy  between  the  ionization  limit  and  the 
level  to  which  the  electron  recombines,  Ejx,  and  the  kinetic  translational  energy  of  the  recombining 
electron.  Etr-  It  is  this  dependence  on  the  electron  translational  energy  that  allows  the  determination 
of  the  electron  temperature  if  the  electrons  have  a  Maxwellian  velocity  distribution.  The  needed 
cross-section  information  describing  this  recombination  is  contained  in  the  so-called  Biberman 
factor./5/  For  this  research,  all  values  for  this  factor  were  taken  from  Hofsaess./6/  At  wavelengths 
above  400  nm,  the  dependence  on  Te  is  very  weak  and  one  can  measure  nc  for  sinelv  ionized 
plasmas  such  as  ours  where  n;  =  ne.  We  chose  to  use  428.5  nm  for  our  measurement.  Our  inferred 
He  was  also  compared  to  that  obtained  from  centerline  corrected/?/  Stark  widths  of  the  H13  line,  and 
agreement  to  10%  was  found.  Alternatively,  with  ne  one  can  make  another  emission  measurement  at 
a  wavelength  below  400  nm  where  the  dependence  on  Te  becomes  strong  enough  to  measure.  We 
chose  270  nm  to  yield  Te  with  an  uncenainty  of  6%.  In  general,  this  method  provides  Te  but  not  the 
gas  temperature,  T^.  However,  for  plasma  conditions  such  as  ours  where  there  are  minimal 
electromagnetic  fields  present,  moderate  radial  temperature  gradients,  and  minimal  radiation  escape, 
ihe  electron  energy  equation  predicts  a  one  kinetic  temperature  plasma  such  that  Ttj  equals  Te- 

5.  Results  and  Discussion 


Emission  measurements  were  made  at  both  the  nozzle  and  test  section  exit  for  a  pure  argon 
plasma  and  then  for  an  argon  plasma  with  varying  amounts  of  either  nitrogen  or  hydrogen  added. 

The  pure  argon  plasma  has  been 
characterized  as  an  LTE  plasma  within 
about  1  cm  of  the  radial  center  at  the  nozzle 
exit  and  as  a  PLTE  plasma  within  about  1.6 
cm  of  the  radial  center  at  the  test  section 
exit./2/  The  diluent  is  added  to  this 
established  experimental  condition.  In  all 
cases,  emission  coefficients  were  measured 
at  715  nm,  430  nm.  428.5  nm.  and  270  nm 
to  yield  argon  nap,  argon  nsp,  total  ng,  and 
Te,  respectively.  Figure  1  shows  rup,  nsp, 
and  ne  for  all  test  conditions  at  the  test 
section  exit.  We  clearly  see  a  drop  in  all 
number  densities  with  the  addition  of  diluent,  although  the  cause  is  not  apparent  from  this  fisure 
alone.  .Vote  that  this  and  several  other  graphs  have  two  abscissas.  The  solid  data  correspond  to  the 
addition  of  nitrogen  into  the  argon  plasma  with  the  per  cent  of  atomic  nitrogen  added  read  on  the 
bottom  abscissa.  The  open  data  correspond  to  the  addition  of  hydrogen  with  the  amount  read  on  the 
lop  abscissa. 


Per  Cant  Atomx:  Hydrogen 


Per  Cent  Atomic  Nitrogen 
Figure  1 :  Absolute  Number  Densities 
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Since  the  undiluted  plasma  is  in  PLTE  at  the  test  section  exit,  the  bound-excited  and  free 
electrons  are  in  Saha  ecjuiiibrium  at  Tg  such  that  Tg  determines  the  ratio  of  nj  to  ng*".  One  can  then 
evaluate  the  ratio  of  nj  to  ng^  for  cases  with  diluent  normalized  to  the  case  without  diluent  to 

investigate  whether  PLTE  is  maintained. 
As  Figure  2  shows,  this  parameter 
decreases  for  increasing  amounts  of 
diluent.  If  PLTE  still  holds,  this  decrease 
implies  an  increase  in  Tg.  However,  as 
will  be  demonstrated  shortly  (Figure  3),  Tg 
acmally  decreases  with  increased  amounts 
of  diluent.  We  therefore  conclude  that  the 
addition  of  diluent  causes  the  plasma  to 
depart  from  PLTE. 

To  quantify  this  departure  from  PLTE, 
Tg  is  required.  .Although  in  theory  the 
ultraviolet  recombination  measurement  can 
directly  measure  Tg,  in  practice  its  uncertainties  -  especially  near  the  radial  center  and  edges  of  the 
plasma  -  preclude  its  use  as  a  sufficiently  precise  technique.  However,  because  the  pure  argon 
plasma  has  been  shown  to  be  in  PLTE  at  the  test  section  exit,  both  the  Boltzmann  temperature 
inferred  from  relative  argon  excited  state  populadons  and  the  line  to  continuum  temperature  inferred 

from  an  excited  state  Saha  equation  equal 
the  electron  temperature.  These  two 
temperature  profiles  agree  well  out  to  about 
1.6  cm  with  both  Tg/S/  and  with  a 
temperature  profile  generated  from  a  fluid 
dynamics  code/8/  constrained  to  agree  at 
the  radial  center.  Further,  the  code's 
output  was  perturbed  by  less  than  2% 
when  the  thermodynamic  propenies  were 
adjusted  to  account  for  the  presence  of  the 
diluent.  Consequently,  the  computer 
generated  shape  is  taken  as  the  electron 
Figure  3:  Entfialpy  and  Electron  Temperatures  temperature  shape  for  all  cases  with  its 

absolute  masnitude  known  for  pure  argon.  To  evaluate  the  magnitude  of  this  profile  for  cases  with 
diluent,  a  calorimetric  energy  balance  was  utilized. 

The  enthalpy  tlux  was  monitored  at  both  the  nozzle  and  test  section  exit.  For  all  experimental 
cases,  the  measured  sensible  fluxes  were  kept  within  3%  of  each  other  at  the  nozzle  exit  through 
adjustment  of  the  plate  input  power.  That  is,  the  input  power  was  adjusted  so  that  the  nozzle  exit 


Per  Cent  Atomic  Hydrogen 


Per  Cent  Atomic  Nitrogen 


Per  Cent  Atomic  Hydrogen 


Rgure  2:  PLTE  Companson  Parameter 
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temperatures  remained  nearly  the  same  in  the  presence  of  the  diluent.  Because  the  8400  K 
temperature  at  the  nozzle  exit  is  sufficient  to  completely  dissociate  the  diluents,  the  corresponding 
total  enthalpy  fluxes  for  the  cases  with  hydrogen  and  particularly  nitrogen  -  because  of  nitrogen's 
ciissociation  enersy  *  are  larger  than  that  for  pure  argon.  Since  the  nozzle  exit  temperatures  are 
then  ail  taken  to  be  the  same,  one  can  infer  Tg  at  the  test  section  exit  for  each  case  with  diluent  by 

comparing  its  measured  enthalpy  drop 
across  the  test  section  including  radiative 
losses/2/  to  the  enthalpy  drop  for  the  pure 
argon  case.  For  example,  at  the  test  section 
exit,  we  calorimetrically  measure  a  sensible 
enthalpy  of  10  kW  for  pure  argon.  For 
every  additional  1  kW  of  calorimetric  loss 
measured  across  the  test  section  resulting 
from  the  addition  of  diluent,  then,  Tg 
would  drop  10%.  This  approach  assumes 
that  there  is  negligible  recombination  of  the 
diluent  between  the  nozzle  exit  and  the  tes: 
section  exit.  An  equilibrium  analysis  of  the 
various  gas  compositions  at  the  inferred  Tg  show  that  this  is  valid  for  all  test  conditions  except  for  the 
largest  addition  of  nitrogen.  Preliminary  modelling  of  this  recombination  process  with  time 
dependency  -  CHEMKIN  -  suggests  that  nitrogen  recombination  would  lag  equilibrium  for  our 
conditions.  Experimental  emission  measurements  of  the  N2'''  1st  negative  and  N2  2nd  positive  band 
lieads  also  show  that  compared  to  equilibrium,  the  emission  is  less  than  expected  indicating  hindered 
recombination  of  atomic  nitrogen.  Plotted  in  Figure  3,  then,  are  the  resulting  test  section  exit 
temperatures  0.5  cm  from  the  radial  center  along  with  the  corresponding  ultraviolet  Tg.  Recalling 
that  Tg  is  6%  uncertain,  generally  good  agreement  is  observed. 

With  Tg  the  deviation  from  LTE  and  PLTE  can  be  quantified.  Under  the  assumption  of  LTE  for 
our  atmospheric  pressure  plasma,  equilibrium  argon  excited  state  populations,  nj*,  and  equilibrium 
electron  number  densities,  ng*,  can  be  calculated.  A  useful  non-dimensional  parameter  which  then 


reflects  the  deviation  from  LTE  is  alpha,  defined  by  the  following: 


a.  4 


“e  = 


* 


V"ey 


(3) 


Plotted  in  fiaure  4  is  alpha  for  the  pure  argon  case  at  the  test  section  exit.  We  see  that  the  plasma 
deviates  by  over  an  order  of  magnitude  from  LTE  for  radial  positions  greater  than  about  1 .3  cm.  but 
that  the  asreement  among  all  the  alphas  implies  PLTE.  This  is  one  of  the  primary  tests  of  PLTE. 
Fisure  5  plots  alpha  for  two  representative  cases  with  diluent,  also  at  the  test  section  exit.  Here  we 
see  more  evidence  of  non-PLTE.  For  both  nitrogen  and  hydrogen,  the  alphas  differ  by  a  factor  of 
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jetween  five  and  ten  demonstrating  the  absence  of  PLTE.  We  also  see  that  the  4p  level  is  closest  to 
LTE  and  the  free  electrons  are  funhest  from  LTE  suggesting  the  diluent's  effect  is  upon  argon’s 
lower  energy  levels. 

To  bener  understand  this  phenomena,  we  need  to  examine  first  why  the  pure  argon  case  deviates 
from  LTE.  Recall  that  we  take  our  emission  measurements  nearly  30  cm  downstream  of  the  center 

of  the  coil  region  where  the  power  is 
deposited.  This  creates  a  recombining 
plasma  as  it  loses  energy  primarily  through 
conduction  and  radiation.  Because  the 
electronic  structure  of  argon  allows  excited 
states  to  collisionally  or  radiatively  return  to 
ground  from  states  that  are  all  over  1 1.5  eV 
in  energy  above  the  ground  state, 
collisional  de-excitation  an  electron  or 
argon  atom  or  ion  is  slow.  Radiative  de¬ 
excitation  to  the  ground  state  is  heavily 
trapped  which  further  hinders  the 
recombination  process.  Diffusion  of 
excited  states  and  free  electrons  from  the  honer  core  to  the  colder  plasma  edges  also  contnbutes  to  the 
deviation  from  LTE  creating  the  largest  deviations  near  the  plasma  edges.  According  to  argon 
collisional-radiative  models,/9-10/  if  ne  is  sufficient,  PLTE  with  the  bound-excited  levels  -  especiallv 
the  higher  energy  states  -  is  expected,  and  this  PLTE  is  indeed  observed  down  to  the  lowest  level 
accessible  through  emission,  the  4p.  Since  eiecffons  recombine  predominantly  with  arson  ions  to 
form  an  excited  state  rather  than  a  ground  state  atom,  for  the  recombining  pure  argon  plasma, 
deviations  from  LTE  are  controlled  by  the  finite  recombination  rates  of  the  argon  excited  states  to 
ground. 

With  the  introduction  of  either  atomic  nitrogen  or  hydrogen,  however,  argon  excited  states  have 
additional  collision  partners  for  transitions  to  ground.  This  type  of  de-excitation.  or  quenching, 
caused  by  molecular  diluents  interacting  with  argon's  4s  level  is  already  well  documentedVl  1-13/ 
Tliese  previous  experiments  suggest  that  the  transitions  in  the  molecular  diluents  which  quench  best 
are  ones  with  an  energy  separation  that  is  less  than  that  for  the  argon  transition  to  ground.  These 
types  of  transitions  result  in  a  relatively  fast  exothermic  two-body  reaction  in  contrast  to  a  slower 
endothermic  one  if  the  diluent  s  transition  has  an  energy  separation  larger  than  argon  s.  Both  atomic 
nitrogen  and  hydrogen  have  at  least  two  such  exothermic  transitions.  Neon,  in  contrast,  has  only 
endothermic  quenching  transitions;  therefore  an  experiment  with  the  addition  of  neon  as  the  diluent 
will  serv  e  as  an  expenmental  check./l  1/  The  following  reactions  are  representative  examples; 

A(4s)  +  H(l)  4-^  A(3p)  +  H(2)  +  1.42  eV  (4) 

A(4s)  +  Nf2s)  A(3p)  -t-  N(2p)  0.70  eV 

A(3d)  -t-  Ne(2p)  -t-  2.37  eV  A(3p)  +  Ne(3s) 
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The  quenching  of  argon  excited  states  by  the  addition  of  certain  diluents  now  explains  the 
observed  reduction  in  a4p  and  asp  as  compared  to  the  pure  argon  case  (Figures  4  and  5).  To 
understand  the  departure  from  PLTE,  finite  electron  recombinadon  rates  and  finite  electron  collision 
rates  must  be  considered.  The  differences  between  Og  and  both  aap  and  asp  suggest  that  the 
electron  recombinadon  rate  is  not  fast  enough  to  follow  the  reladvely  quick  quenching  of  argon 
excited  states.  The  differences  between  oup  and  asp  suggest  that  the  electron  collision  rate  with  the 

argon  excited  states  is  also  not  fast  enough 
-  again  relative  to  the  diluent's  quenching 
rate  -  to  maintain  even  a  PLTE  among  just 
the  excited  states.  This  model  implies  a 
competition  between  electron  collisions 
with  argon  excited  states  attempting  to 
maintain  PLTE  with  the  over-populated 
electrons,  and  diluent  collisions  attempting 
to  quench  these  excited  states  in  the 
direction  of  approximate  equilibrium  with 
the  ground  state.  For  analytical  purposes, 
we  define  nj*c  to  be  a  given  argon  excited 
state’s  population  in  equilibrium  with  the 


Per  Cent  Atomic  Nitrogen 
Rgure  6;  Quenching  Parameter 


free  electrons  at  Te  and  nj*h  to  be  a  given  argon  excited  state’s  population  in  equilibrium  with  the 
ground  state  at  Te-  Plotted  in  Figure  6  is  a  parameter  which  combines  these  defined  populations  with 
tho.se  experimentally  determined  to  quantify  the  degree  of  quenching.  This  parameter,  however, 
better  describes  deviations  from  PLTE  and  relative  degree  of  quenching  than  it  descnbes  absolute 
quenching.  Alpha  is  best  for  determining  whether  equilibrium  with  the  ground  state  exists.  In  any 
case.  Figure  6  clearly  shows  that  the  degree  of  quenching  at  the  test  section  exit  is  increased  for 
increasing  amounts  of  diluent,  and  that  even  0.07%  atomic  hydrogen  causes  deviations  from  PLTE. 
We  also  see  that  as  expected  the  4p  level  which  is  closer  in  energy  to  either  of  the  possibly  quenched 
argon  levels  is  more  strongly  depopulated  than  is  the  5p  level.  Also  plotted  is  one  data  pair  for  5% 
neon  which  shows  minimal  quenching  thereby  supporting  the  idea  that  the  diluent  should  have  a 
exothermic  transition  with  respect  to  the  relevant  transition  in  argon. 


6.  Conglusi.Qn 


Experimental  measurements  of  argon  e.xcited  state  densities  and  electron  number  densities 
indicate  that  a  PLTE  exists  for  a  recombining  pure  argon  plasma  such  that  the  free  and  bound-excited 
electrons  are  in  equilibrium  though  over-populated  with  respect  to  the  ground  state.  The  addition  of 
either  niffogen  or  hydrogen,  however,  in  amounts  as  small  as  0.1%  cause  deviations  from  this 
PLTE.  Both  diluents  have  an  exothermic  transition  with  an  excited  state  in  argon.  Experiments  with 
neon  which  has  an  endothermic  coupling,  in  contrast,  demonstrated  negligible  quenching. 
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